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ABSTRACT 


This is the third in a series of reports on the use of tracking information 
from NASA's Syncom n satellite to define the longitude-dependent gravity 
field of the earth. This set of orbit data shows that the westward drift of 
Syncomll decelerated from 0.480 degree s/day at a nodal longitude of -174.35° 
in early July 1964, to 0.446 degrees/day at a longitude of 159.53° in late Aug- 
ust 1964. Further westward drift of the satellite resulted in an acceleration 
of the westward drift rate to 0.501 degrees/day when Syncom n reached a 
mean longitude of 134.22° in early November 1964. On the assumption that 
this four-month drift of Syncom II was sensitive to only second-order 
longitude-dependent earth gravity (associated with the ellipticity of the 
earth’s equator), as well as sun, moon and ordinary earth zonal gravity 
forces, the magnitude and orientation of such equatorial ellipticity is calcu- 
lated as 


J 22 = - (1>86 ± . 14) x icr 6 , 

which corresponds to a difference in major and minor equatorial radii of the 
earth of 


a Q - b Q » 234 ± 18 ft , 


and 


\ 22 = - (16.6 + 2.9)° , 

which locates the longitude position of the major equatorial axis with respect 
to Greenwich. 


ii 
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THE EQUATORIAL ELLIPTICITY OF THE EARTH 
AS SEEN FROM FOUR MONTHS OF SYNCOM II DRIFT 
OVER THE WESTERN PACIFIC 


by 

C. A. Wagner 

Goddard Space Flight Center 


INTRODUCTION 

This report is one of a series on the use of the tracking data from Syncom II, the first opera- 
tional 24-hour satellite, to define the exact nature of the earth's longitude -dependent gravity field. 
Nineteen orbits for this satellite, calculated from range and range rate observations at the God- 
dard Space Flight Center, have been used to fix the longitude location for Syncom n during a four- 
month drift period in the summer and fall of 1964. As in the previous GSFC orbit determinations 
for Syncom n, longitude -dependent earth gravity was ignored in the calculation of these orbits 
and a new determination of elements was required about every week to maintain sufficiently accurate 
position data for the satellite. 

Previous reports in this series have covered two three-month "free" drifts of this satellite 
over Brazil between 55° West and 63° West (Reference 1), and a two-month free drift period over 
the Central Pacific between 117° West and 163° West (Reference 2). During these long periods 
between gross orbit change maneuvers initiated by onboard propulsion, this nearly synchronous 
satellite has remained for sufficient time over limited longitude sectors of the earth to accumulate 
noticeable drift perturbations caused by small longitude mass anomalies within the earth. 

The first gravity experiment with this satellite (Reference 1) tested the sensitivity of its ob- 
served accelerated drift as a function of time in a limited longitude region with respect to the 
simplest kind of longitude-dependent gravity (that associated with an earth whose equator is ellip- 
tical, not circular). The second gravity experiment with Syncom n (Reference 2) tested the sensi- 
tivity of its observed change of longitude drift rate as a function of mean longitude location over a 
much longer arc with respect to the same kind of earth gravity. 

With minor differences (which could be rationalized as being due to the "noise" of these ex- 
periments) these two independent analyses of widely separated Syncom n drift arcs came to the 
same conclusion— that the earth's longitude gravity field has at least one component (dominant on 
24-hour satellites) represented by an elliptical equator whose magnitude gives a second order 
tesseral constant near 


j 22 = - 1.7 x icr 6 , 
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which represents a difference in major and minor axes of about 65 meters or 213 feet, and a major 
axis orientation near 


The close agreement of the two previous experimental results appeared to confirm the basic as- 
sumption of these analyses, that longitude earth gravity components of higher spacial order than 
that represented by equatorial ellipticity are of negligible long term influence on the high altitude 
24 -hour satellite. 

This report examines the observed drift rate changes of Syncom n over a free drift arc be- 
tween 171° West and 133° East. This new drift arc was initiated by the pulsing of onboard cold 
gas jets on 4 July 1964 which resulted in a reduction of Syncom ITs westward drift rate from 0.7 
degrees/day (at the end of the arc analyzed in Reference 2) to about 0.5 degrees/day. 

The analysis of this new arc data for longitude -dependent earth gravity utilizes exactly the 
same method as developed in Reference 2 and proceeds upon the same basic assumptions as in 
both previous analyses, namely: 

1. Only second-order longitude-dependent earth gravity (associated with an earth whose 
equator is elliptical) was sensed over the long term drift of Syncom n, in this case from 
July to November 1964 

2. All perturbations on Syncom n not due to gravity during this period were negligible 

3. The earth zonal, sun and moon gravity effects on Syncom n drift can be adequately repre- 
sented over the span of a few months by a small bias in the reduced results from the 
actual data, which can be estimated from simulated Syncom n trajectories numerically 
integrated. 

The reader may wish to consult References 1 and 2 for further background on the source of the 
method of analysis used here. 


THEORY OF DATA REDUCTION 

The square of the net drift rate of a 24-hour satellite in a second-order longitude gravity 
field is given (Equation 3 in Reference 2) as 

(\) 2 = - A 22 cos 2(\ -\ 22 ) + c, , (1) 

where C x is a constant determinable from initial or other known or calculable drift conditions. 

The longitude location of the 24-hour satellite is \ 22 is the longitude location of the major axis 
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of the elliptical equator, and A 22 is an earth- gravity -orbit constant given by 


= - 72t7 2 J\ 


+ 1 


rad 

(sid day) 2 


( 2 ) 


In Equation 2, J 22 is the second-order longitude gravity constant in the earth's spherical harmonic 
potential field (Appendix A), R 0 is the mean equatorial radius of the earth and a s is the semimajor 
axis of the nearly circular orbit synchronous satellite of orbit inclination i s . In Reference 3 it is 
shown that Equation 1 holds strictly for a synchronous period, circular orbit satellite in an earth 
field without higher order longitude gravity. For the inclined circular orbit satellite, k is inter- 
preted as the geographic longitude of the equator crossing. The drift rate of this crossing longi- 
tude is k. As long as k is sufficiently small (e.g., k < 2 degrees/day) and the eccentricity of the 
orbit is not excessive (e.g., e < .01), Equation 1 is still a very good approximation of the law of 
drift motion for the 24-hour satellite. For somewhat more restrictive mean drift rates and ec- 
centricities than those stated above, the adequacy of Equation 1 has been confirmed by many nu- 
merically calculated 24-hour partial drifts as well as by the second integral of the pendulum drift 
motion which stems from Equation 1. The simple order of magnitude arguments of Reference 3 
are given in support of the stated probable upper limits on k and e, for which Equation 1 holds 
with reasonable accuracy. 

For the analysis of the "free" Syncom II drift in July-November 1964 it is merely noted that 
k < 0.51 degrees/day and e < 0.001 applies throughout. To minimize the small error involved in 
assuming A 22 is a constant in the drift, it is assumed that A 22 is determined from average values 
of the relevant orbital parameters (i s and a s ) over the complete drift period. The arbitrary con- 
stant of the motion, c i? will be calculated in the analysis as a quasi-mean, constant drift rate for 
the period about which the longitude -gravity -caused perturbation in the drift rate occurs. Equa- 
tion 1 can be expanded and rewritten as 


CM 2 = C, + C 22 Fj (A) + S 22 F 2 (A) 


where 


(3) 


F, (M 


367 T 2 



(cos 2 i s + 1 ) cos 2 k 


F 2 (A) 


3677 2 



2 

I (cos 2 i s + 1 ) sin 2k 


C 22 J 22 cos 2 ^22 


and 


(4) 

(5) 

( 6 ) 


S 22 - J 22 sin2X 22 


(?) 
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Assume three or more sets of X, X values over this free drift arc. Also assume average values 
for the slightly varying a s and i s parameters of the orbit. Then, "best" values of the constants 
C 2 , C 22 , and S 22 can be found which minimize, for example, the sums of the squares of the devia- 
tions of (X ) 2 (observed) from (X ) 2 as calculated from Equation 3 with this "best set" of drift par- 
ameters, Cj, C 22 , and S 22 . This rational data- smoothing procedure is known as "least squares" 
fitting. 

Given the least squares values of c 22 and § 22 , from Equations 6 and 7 we can calculate "best 
values" of j 22 and X 22 from these derived ellipticity parameters as 


and 


1 22 


22 


1 

n> 

5* 

+ s 2 2 2 ) 

i , 

/-s 22 

2 tan 


-r 

\ C 2 2 


1/2 


(8) 

( 9 ) 


The negative signs in numerator and denominator of Equation 9 appear because, from the definition 
of k 22 as the major equatorial axis and the form in which the earth potential function is written in 
this analysis (Appendix A), j 22 must be negative. Thus these mandatory negative signs cannot be 
arbitrarily cancelled if the quadrant of X 22 is to be properly placed. 


UNSTABLE EQUILIBRIUM LONGITUDE 

In Figure 1 it is observed that Syncom n drift in July-November 1964 passes across a 
longitude X e where 


dX 
dX e 


0 . 


Since 


dX dX dt X 
dX “ dt dX ” \ 

and X e was finite, this implies X e = 0, or X e is an equilibrium longitude in the drift regime of this 
24-hour satellite. Since the drift rate slowed approaching X e (in late August 1964) and increased 
afterwards, X e is evidently a longitude of unstable equilibrium. As such, for the second- 
order longitude-dependent gravity field assumed in this analysis, it is shown in Reference 1 
that x e (unstable) is just one of the longitude extensions of the major axis (x 22 ) of the ellip- 
tical equator. 

Thus, this unstable equilibrium longitude (one of four equilibrium longitudes around the equa- 
tor and an important parameter for the design of operational synchronous satellite systems) 
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Figure 1 — Drift rate of Syncom II (summer and fall, 1964) between 188° and 134° East. 


was found directly from X 22 as measured in this Syncom n - Central Pacific gravity experiment. 


ANALYSIS OF NINETEEN SYNCOM II ORBITS IN THE SUMMER AND FALL 
OF 1964 TO DETERMINE THE ELLIPTICITY OF THE EARTH’S EQUATOR 

Table B-l in Appendix B lists nineteen T, orbits M for Syncom II, about a week apart, which were 
calculated from returned tracking information on the satellite at Goddard Space Flight Center. 
These orbits run from 7 July 1964 (when the mean westward drift rate was slowed from 0.7 degrees/ 
day to 0.5 degrees/day at about 171° West) to 2 November 1964, when the equator crossing of the 
satellite was at about 133° East.* We can obtain an estimate of the longitude drift rate of the 
satellite between orbit epochs by simply talcing the ratio of the difference in ascending equator 
crossing longitudes to the difference in time between these near epoch crossings. If this mean 


'Table 1 presents a list of ascending equator crossings for Syncom II nearest the epochs of the GSFC-calculated elements in Ap- 
pendix B. These have been derived by numerical trajectory generation from these elements using the earth -gravity, sun and moon 
constants identical to those used in the GSFC orbit determination program. 
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drift rate between crossings is then assigned to the mean longitude in that interval, a small longi- 
tude error will be made which appears to be negligible. If the drift interval is less than about 
10° or the time interval less than about 10 days, Appendix C of Reference 2 shows that the error 
made in this estimate should be negligible. In addition, the efficacy of this mean drift and drift 
rate technique is well proven by the results of the data reduction on the simulated Syncom II tra- 
jectory in Appendix C of this report. (In the next section will be found another technique for drift rate 
estimation which may prove useful in future studies of24-hour satellite motions.) Table lalso pre- 
sents this simply derived drift rate as well as inclination and semimajor axis data at ascending equa- 
tor crossings for Syncom n in July-November 1964. At the bottom of Table 1 arefound average values 
of i s and a s in this drift period. An "observation" on the drift rate squared as a function of longitude 
(according to the theory of Equation 1 or Equation 3) is then taken by the following set of data: 

1. the average i s for the entire drift period, 

2. the average a s for the entire drift period, 

3. a set of mean (bracketed) X and X values. 


With this data Fj ( x ) and F 2 ( X ) can be calculated from Equations 4 and 5 and a condition 
equation on the "observed" (\) 2 can be written to satisfy Equation 3. We then assume that each 
such condition equation on the observed data also contains a small additive and varying amount of 
error e (the so-called residual for that condition equation). The method of least squares estima- 
tion of the drift and ellipticity parameters C 19 C 22 , and s 22 finds those parameters (c x , C 22 , S 22 ) 
which minimize the sums of the squares of the e f s. These residuals are determined by the dif- 


ferences of (>0 o 2 bserved and C x + C 22 F, (^) observed 


+ S 22 F 2 (\) 


obs e rved 


for the complete set of obser- 


vations. One then can assume for convenience that the residuals are random and normally distri- 
buted among the observations with a constant variance a 2 and a mean value zero. In this case it 
can be shown (Reference 4) that c l? C 22 , and § 22 are also random and normally distributed with 
variances proportional to cr £ andeasily calculated from the so called’ 'normal equations" of condition 
(see also References 2 and 5). The (bracketed) drift and drift rate observations of Table 1 have 
been inputed into a revised IBM Share program (Reference 6) which sets up and then solves the 
"normal equations” of condition corresponding to Equation 3. 


As outputs in this program (checked for accuracy on the data in Reference 2 which was pre- 
viously worked up on a desk calculator), estimates of cr e (called s e ), <t(c 22 ) (called s(c 22 )), and 
cr(s 22 ) (called s(s 22 )), as well as C %9 C 22 , and S 22 have been obtained and are listed at the bottom 
of Table 1. 


The results of this statistical analysis of the drift data on Syncom n (according to the theory 
of Equation 1 or Equation 3) during July-November 1964 can be summarized by the following least 
square estimates and standard deviations: 


- ( 1. 6123 ±0. 1439) x 10' 6 

(10) 

(1.0683 +0.1140)x 10‘ 6 . 

(11) 
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Table 1 


Syncom n Ascending Equator Crossings and Related Drift Data, Summer and Fall 1964.* 
Initial Orbit Epoch (yr-mo-day-hr, universal time): 64-07-04-2.0 


Crossing 

Number 

(J) 

Date, 

1964 

(mo- day) 

Time, Tj 
(days from 
Jan. 0.0, 1964) 

\T, 

((Jays) 1 

longitude,* Lj 
(degrees) 

\L, 

L f L |-l 

(degrees) 

Observed 
Longitude 
Drift Rate,* 
M/AT 
(Jegrees/day) 

(ncfinatron 
at Crossing, 
l s 

(32.0°Mtata) 

Semimajor 
Axis af 
Crossing, a s 
(6.616 + 
data x 10,' 4 
earth radii) 

Square of 
Observed 
Longitude 
On ft Rate, 

Y (observed) 

{ fad/sid dayf 

Square of 
Computed 
Longitude 
Drift Rate,** 
Y (computed) 
(latl'sid day) 2 

Square of Residual 
Longitude Drift 
Rale, 

Yfobs.) -Y(comp.) 
(rad/sid day) 2 

1 


186.5598 


-171.249 



.541 

4.999 





7-6 



(-171.96) 


(-.47435) 



6.817x10-3 

7.056x 10-s 

-2.396x 10-6 

2 


189.5555 

2.9957 

-172.670 

-1.421 


.537 

8.957 





7-11 



(-174.35) 


(-.4804) 



6.992 

6.893 

9.886xl0' 7 

3 


196.5455 

6.9900 

-176.028 

-3.358 


.500 

7.147 





7-18 



(-177.90) 


(-.4696) 



6.681 

6.675 

6.146x 10- 8 

4 


204.53365 

7.98815 

-179.7795 

-3.7515 


.516 

5.120 





7-25 



(178.82) 


(-.4678) 



6.630 

6.502 

1.279 x 10" 6 

5 


210.5249 

5.99125 

177.418 

-2.8025 


.500 

6.072 





8-1 



(175.80) 


(-.4628) 



6.489 

6.369 

1. 196x 10' 6 

6 


217.51435 

6.98945 

174.183 

-3.235 


.476 

6.558 





8-8 



(172.595) 


(-.4544) 



6.255 

6.258 

-2.348x 10-0 

7 


224.5039 

6.98955 

171.007 

-3.176 


.444 

2.387 





8-15 



(169.44) 


(-.4498) 



6.129 

6.179 

-4.950xl0- 7 

8 


231.4932 

6.9893 

167.8635 

-3.1435 


.409 

4.756 





8-22 



(166.28) 


(-.4520) 



6.1895 

6.132 

5.791 x 10- 7 

9 


238.4828 

6.9896 

164.704 

-3.1595 


.443 

4.324 





8-29 



(163.14) 


(-.4467) 



6.045 

6.117 

-7. 134 x 10~ 7 

10 


245.4718 

6.9890 

161.582 

-3.122 

1 

1 

.397 

2.705 





9-6 



(159.53) 


(-.44635) 



6.036 

6.139 

-1.03 1 x 10-6 

11 


254.4580 , 

8.9862 

157.521 

-4.011 


.372 

0.392 





9-13 


! 

(156.20) 


(-.4564) 

! 


6.311 ! 

6.197 

1. 140xl0' 6 

12 


260.4489 ' 

5.9909 

154.8365 

-2.7345 


.360 

4.606 ' 





9-19 j 



(153.47) j 


(-.4571) 


I 

6.330 

6.270 

5.984 x 10‘ 7 

13 


266.4398 ! 

5.9909 

152.098 

-2.7385 


.302 

3.414 





9-26 



(150.49) 


(-.4608) 



6.433 

6.376 

5.645x 10- 7 

14 


273.4294 

6.9896 

148.877 

-3.221 


.327 

6.017 





10-3 



(147.29) 


(-.4534) 



6.228 

6.5195 

-2.916x 10‘ 6 

15 


280.4188 

6.9894 

145.708 

-3.169 


.321 

4.898 





10-10 



(144.055) 


(-.4730) 



6.778 

6,693 

8.513 x 10 ‘ 7 

16 


287.4085 

6.9897 

142.402 

-3.306 


.312 

7.049 





10-17 i 

i 


(140.51) 


(-.4749) 


1 

6.830 

6.913 

-8.370x 10- 7 

17 


295.3969 

7.9884 

138.609 

-3.793 


.290 

3.824 





10-24 



(137.17) 


(-.48155) 



7.025 

7.147 

-1.222 x 10'6 

18 


301.3880 

5.9911 

135.724 

-2.885 


.250 

10.382 





10-30 



(134.22) 


(-.5012) 



7.610 

7.373 

2.376x IQ ' 6 

19 


307.3798 

5.9918 

132.721 

-3.003 


.249 

5.259 





Average < s = 32.397°; Average « s = 6.6165203 


*The crossing longitudes without brackets are the first ascending equator crossings of Syncom II after the epoch's of Table B-l as 
reported by the Goddard Space Flight Center Tracking and Data Systems Directorate. The longitude drift rates within brackets 
are the mean longitude rates between successive epoch crossings assumed to occur at a longitude midway between these 
crossings. 

♦♦Computed from Equation 3 using the values of C, } C 22 and S 22 below. 

Results from the "Least Squares" estimation of the parameters C 1 , C 22 and S 22 in Equation 3 for this actual Syncom II drift: 

C, = 8.7863x10-5 rad/sid. day 2 S(C 22 ) = 1.439xl0‘ 7 

C 22 =-1.6123xl0- s S(S 22 )= 1. 140x 10' 7 

S 22 = 1.0683 x IQ -6 S f (standard error of estimate) = 1.437x 10 6 rad/sid. day 2 



~ _ rad 

C, = 8.7863 x 10" 5 — ; 

1 (sid day) 2 


1.437 x 10' 6 


rad 
(sid day) : 


(12) 

(13) 


Figure 1 presents a graph of the "observed” longitude versus longitude drift rate data for 
Syncom n in this period together with a theoretical curve calculated from Equation 3 with the above 
estimates C x , C 22 and S 22 . In addition, a dashed, broken line gives the results of a numerically 
calculated (simulated) Syncom n trajectory identical to that in Appendix C except that longitude- 
dependent earth gravity was not included in the calculation (see Appendix C, Table C-2). The dis- 
crepancy between this almost level line and the actual "bowed" data is explained with a fairly 
reasonable tolerance by the earth’s longitude gravity field as represented through Equation 3 by 
the C 22 and S 22 parameters in Equations 10 and 11. 

The second-order earth longitude gravity constants J 22 and \ 22 (Appendix A) corresponding 
to Equations 10 and 11 are calculated from Equations 8 and 9 respectively as: 


J 22 = - ( 1. 93 ± 0. 1 3) x 1CT 6 


(14) 


and 


£ 22 = ~ ( 16.8 ± 2 . 8 )° . ( 15 ) 

(See Appendix E of Reference 2 for the treatment of the estimation of the standard deviations of 
j 22 and \ 22 from c 22 , S 22 and their standard deviations.) 

As the dashed line in Figure 1 shows, the sun, moon and earth zonal gravity fields alone have 
some small variable influence on the drift rate of Syncom II in this four month period. Therefore, 
we treat the sensed earth gravity constants in Equations 14 and 15 as unadjusted for these 
effects (as well as second-order model error due to the approximate nature of Equation 1) and as- 
sume they contain a small bias such that: 


J 22 (adjusted for sun, moon, - J 22 (unadjusted) + J 22 (bias) (16) 

earth zonal and second 
order longitude gravity 
error) 


\ 22 (similarly adjusted) - \ 22 (unadjusted) + k 22 (bias) 


(17) 


It is noted that we do not make any attempt to adjust over this limited amount of rather noisy data 
for the presence of an unknown amount of higher order earth longitude gravity (refer to the next 
section). The simulated Syncom n trajectory in Appendix C with crossing data reasonably close in 
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time and longitude to the actual Syncom n data provides a means both to test the simple drift rate 
theory of Equation 1 and to arrive at an estimation of the biases in Equations 16 and 17. 

From this simulated Syncom n trajectory (July -November 1964) reduced by the same methods 
as were used to reduce the actual trajectory (with the sun, moon and other influences treated as 
noise), it is estimated (Appendix C) that the numerically calculated drift was compatible with appar- 
ent earth longitude gravity constants J 22 and a. 22 , where 

J 22 (apparent) = - (1.75±.04)x 10‘ 6 (18) 

and 

\ 22 (apparent) = - (18.2±.6)° . (19) 

Since J 22 = - 1.68 x 10 6 and ^ 22 = - 18.0° were the actual constants used in the trajectory programs 
that generated the crossing data, we observe: 

1. that accumulated sun, moon, earth zonal and second-order model longitude gravity error, 
on the theory of Equations 1 or 3, is small over the four-month drift period; and 

2. that an estimate of the bias of such accumulated effects is contained (from Equations 18 
and 19) in Equations 16 and 17 with 

J 22 (adjusted) = - 1.68 x 10~ 6 and J 22 (unadjusted) = - 1.75 x 10" 6 , 

\ 22 (adjusted) = - 18.0° and k 22 (unadjusted) = - 18.2° f implying that 
J 22 (bias) = - (1.75 ±0.04)x 1(T 6 - (- 1.68x 10" 6 ) 

= - (0.07 ±0.04) x 10" 6 , (20) 

£ 22 (bias) = - (18.2+0.6)° - (-18.0°) 

- - ( 0. 2 ± 0 . 6 )° . ( 21 ) 

Applying these calculated biases to Equations 14 and 15, we arrive at adjusted values of the ac- 
tually measured Syncom n gravity constants for the July-November 1964 drift: 

J 22 (measured and adjusted = - (1.86 ±0. 14) x 10“ 6 (22) 

for sun, moon and known 
earth model effects in 
Syncom II western Pacific 
drift) 

k 22 (measured and similarly = - (16.6 ±2.9)° . (23) 

adjusted) 
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The value in Equation 22 corresponds to a difference in major and minor equatorial axes of the 
geoid of 


a Q - b Q = 233 ± 17 feet . (24) 

The value in Equation 23 senses the longitude location of the major equatorial axis of the elliptical 
equator at (16.6 ± 2.9)° West. 

The discussion covering the theory of data reduction shows that the other extension of the 
equatorial major axis at (163.4 ± 2.9)°, representing a dynamically unstable equilibrium longitude 
for Syncom n, can be considered to have been directly measured in this gravity experiment. If 
higher order earth longitude gravity effects are reasonably small as they appear to be (refer to 
the next and the last sections), this unstable equilibrium longitude measurement will be valid for 
all 24-hour, near-circular orbit satellites if small epoch-variable sun and moon effects are ignored. 

DISCUSSION 

The most gratifying aspect of the results of this straightforward analysis of nineteen new 
consecutive Syncom n updated orbits is, once again, the close agreement with the previous gravity 
experiments on this satellite over Brazil and the Central Pacific. We have now looked at 54 orbits 
for Syncom n between 133° East and 305° East longitude for sensitivity to second order earth tes- 
serai gravity. Between 297° East and 305° East the experiment measured longitude gravity con- 
stants (Reference 1) 


J 22 = - ( 1. 70 ±0.05) x 10“ 6 


(25) 


and 


\ 22 = - (19.0 ± 6) ° . 


(26) 


Between 245° East and 197° East the Syncom n gravity experiment measured (Reference 2) 


J 22 = - (1.72 ±0.22) x 10' 6 (27) 

and 

\ 22 = - (17.1 ±4.8)° • (28) 

In this new experiment between 133° East and 189° East the measured constants are 


J 22 = - (1.86 ±0.14) x 10~ 6 (29) 

and 

\ 22 = - (16.6 + 2.9)° . (30) 
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It appears, from a comparison of these numbers with their standard deviations, that one of 
two conclusions can be supported from the Syncom n study of second order longitude gravity thus 
far. They are: 

1. longitude earth gravity higher than second order has negligible influence on the 24-hour 

satellite (the second-order gravity constants are roughly J 22 = - (1.8 db .15) x 10" 6 and 

X 22 = - (18 ± 4) °), 
or 

2. longitude earth gravity higher than second order has a small influence on the 24-hour 
satellite which accumulates about a 10 percent increase in the maximum longitude force 
over the western Pacific as compared to the maximum force over western Brazilian 
longitudes. 

In view of the study in Reference 7 (see also Figure A1 in Reference 1) on the probable influence of 
higher order gravity on the Syncom n Brazil measurement, the author is inclined towards the 
second conclusion above. In that study it is recalled that gravity data from lower altitudes showed 
unambiguously that the Syncom n Brazil measurement of J 22 was likely to have been underestimated 
by as much as 25 percent. In progress is a comprehensive analysis of all Syncom II and Syncom m 
drift data thus far which should clarify this situation. 

In the study of Syncom II Pacific drift so far, longitude rates have been estimated from dif- 
ferenced crossing data between independently determined "orbits." An alternate approach to the 
estimation of sets k and k during the drift of Syncom II has been tried and rejected for this study 
but may find use as a supplement to the mean drift method in future studies of the drift of Syn- 
com 13 and Syncom Ht. This method is to use the drift rate shown by the first few observed as- 
cending equator crossings for the satellite after each orbit epoch. In theory the use of such drift 
rate data would almost double the number of independent observations analyzed in each free drift 
arc. The approach, thus far, has been to use only position data from each orbit. But an orbit is 
defined (in theory) by independently given velocity as well as position data. Using the first few 
equator crossings after each epoch as the net drift rate observation for that orbit is equivalent to 
using the velocity data in that orbit determination. In turn, the accuracy of this approach depends 
primarily on the accuracy of the semimajor axis determination for Syncom U. Unfortunately, 
this determination appears to be often too poor to provide sufficient precision to the measurement 
of drift accelerations in this experiment (see, for example, Figure 1 of Reference 2). Neverthe- 
less, the use of such "short arc" data to supplement the "smoothed" rate data derived from the 
longer arcs between "orbits" appears promising over a good portion of the Syncom II arcs already 
studied. 

CONCLUSIONS 

Study of the "free" longitude drift acceleration of Syncom II during July-November 1964 over 
the western Pacific showed sensitivity to the following parameters of earth equatorial ellipticity: 

J 22 = - (1.86 ±0.14) x IQ’ 6 (31) 
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(corresponding to a difference in major and minor equatorial axes of: a 0 -b Q = 233 ± 17 ft), and 


X 22 = - (16.6 ±2.9)° (32) 

(locating the geographic longitude of the major equatorial axis). 

The above results when compared with previous gravity drift studies of Syncom II appear to 
be the first measured indication of the small influence of higher order earth longitude gravity on 
the 24-hour satellite. 

Syncom n in its drift over the western Pacific passed over a dynamically unstable equilibrium 
longitude in the geographic regime of the 24-hour satellite located at 

\ e - 163.4 ± 2.9° . (33) 

These independent experimental results again confirm the conjecture in the first of this series 
of Syncom n gravity investigations that longitude earth gravity of higher than second order has a 
small influence on the long term drift of the high altitude 24-hour satellite. 

A comprehensive analysis of all Syncom EE and EH drift data through 1964 is now underway by 
the author to define the earth T s longitude gravity field at 24-hour altitudes for longitudes from 65° 
to 305° East. 


(Manuscript received June 3, 1965) 
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Appendix A 


Earth Gravity Potential and Force Field Used in This Report 

The gravity potential used as the basis for the data reduction in this study is the exterior po- 
tential of the earth derived in Reference A-l for geocentric spherical coordinates referenced to the 
earth T s spin axis and its center of mass. The infinite series of spherical harmonics is truncated 
after j 44 . The nontesseral harmonic constants j 20 , j 30 , and j 40 are derived from Reference A-2. 

The earth radius R 0 used in this study is 


R 0 = 6378.388 km . 

The earth's gaussian gravity constant used is 

M E = 3.9862677 x 10 s km 3 /sec 2 . 

Neither of these values, taken from Reference A-3, nor the "zonal geoid" of Reference A-2 is felt 
to be the most accurate known to date. They are the values used by the GSFC Tracking and Data 
Systems Directorate to calculate the orbit elements of Syncom n from radar and Minitrack observa- 
tions. They were chosen to insure consistency between the data of this study and these published 
orbits, inasmuch as the "triaxial" reduction for which this study has been undertaken is not sig- 
nificantly affected by the probable errors in these values. The second-order tessera! harmonic 
constants used in the simulation studies were 


J 22 = - 1.68 X 10“ 6 , 

^22 n ” 18 degrees . 

These are the values shown (or implied) on the "tesseral geoid" below (for the J 22 harmonic). At a 
later point in the analysis, the slightly different values reported in the abstract were estimated. 

The most accurate "zonal geoid" is probably that of Kozai (Reference A-4), with the following earth 
constants: 

R 0 = 6378.2 km , 

J 22 ° 1082.48 x 10" 6 , J 30 = ~ 2.56 x 10“ 6 , J 40 = - 1.84 x 10" 6 , 

M e = 3.98603 x 10 5 km 3 /sec 2 . 
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The earth’s gravity potential (to fourth order, probably sufficient to account for all significant 
longitude perturbations on a 24-hour satellite) is given by Equation A-l, which may be illustrated 
(following Reference A-4, Appendix B, with the zonal constants of Reference A-2) as follows: 


NORTH 
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or 


F r - 7l{-l + (R 0 / r ) 2 [ 3 / 2 J 20 (3»in»*-l) + 9J aa co. a *co«2(\-\„) 

+ 2 ( R o/ r ) J 30 (5 sin 2 0- 3 ) (sin0) + 6(R 0 /r) J 31 (S sin 2 0- l) cos 0 cos (A. ~^ 31 ) 

+ 60(R 0 /r) J 32 cos 2 0 sin0 cos 2[K-K 32 ) +60(R Q /r) J 33 cos 3 0 cos 3(^-\ 33 ) 

+ 5/8(R 0 /r) 2 J 40 (35 sin 4 0 - 3Osin 2 0 + 3 ) 

+ 25/2(R 0 /r) 2 J 4 , (7 sin 2 0 - 3 ) cos 0 sin 0 cos (X. - \ 4J ^ 

+ 75/2(R 0 /r) 2 J 42 (7 sin 2 0 - l) cos 2 0cos 2(\ - \ 42 ) 

+ 525(R 0 /r) 2 J 43 cos 3 0 sin 0 cos 3(\-\„ 3 ) + 525(R„/r ) 2 J 44 cos 4 0 cos 4(\-\ 44 )] j- • (A-3) 

F A = 71 ( R o / r ) 2 { 6J 22 cos 0 sin 2(\ -\ 22 ) + 3/2(R 0 /r) J 3 , (s sin 2 0~ l) sin (\-\ 31 ) 

+ 30(R 0 /r) J 32 cos 0 sin 0 sin 2(\ -\ J2 ) +45(R 0 /r) J 33 cos 2 0 sin 3(\ -\ 33 ) 

+ S/2(R 0 /r) 2 J 41 (7 sin 2 0 - 3 ) sin 0 sin (\ - \ 4 , ) + 15(R 0 /r) 2 J 42 (7 sin 2 0- l) cos 0 sin 2(* - \„ 2 ) 
+ 315 (R 0 / r) 2 J 43 cos 2 0 sin 0 sin 3^ - A- 43 ) 

+ 4 2O(R o / r ) 2 J 44 co s 3 0sin4(\-\ 44 )} • (A-4) 

^ e r 

F 0 “ T 2 ^ ( R o/ r ) 2 l 3 J 2 o sin ^ cos <£ + 6J 22 cos <£ sin 0c°s 2(x ~ X 22 ) 

- 3/2(R 0 /r) J 3Q (5 sin 2 0 - l) cos 0 + 3/2(R 0 /r) J 31 (l5sin 2 0- ll) sin 0 cos ( \ - K n ) 

+ 15 ( R o/ r ) J 32 ( 3 sin2 1 ) cos <P cos 2(\ ^ X 32 ) 

+ 45 ( R o/ r ) J 33 c °s 2 0 sir\<fi cos 3(\-\ 33 ) - 5/2(R Q /r ) 2 J 4Q (7 sin 2 <£- 3 ) sin0cos0 
+ 5/2(R 0 /r) 2 J 41 (28 sin 4 0- 27 sin 2 0+3) cos ( A. - \ 4J ) 

+ 30(R 0 /r) 2 J 42 (7 sin 2 0 - 4 ) cos0sin0cos 2(\-X 42 ) 

+ 105(R 0/ /r) 2 J 43 (4 sin 2 0- l) cos 2 0 cos 3(\-\ 43 ) 

+ 42°(R 0 / r ) 2 J 44 co s 3 0sin0c os 4(X-A. 44 )} * (A-5) 
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The actual sea level surface of the earth is to be conceptualized through Equation A-l as a sphere 
of radius 6378 km, around which are superimposed the sum of the separate spherical harmonic 
deviations illustrated. To these static gravity deviations, of course, must be added a centrifugal 
earth rotation potential deviation at the earth T s surface, to get the true sea level surface (Reference 
A-l). 
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Appendix B 

Orbital Elements for Syncom II During the Summer and Fall of 1964 


Table B-l gives the elements for Syncom II during July-November 1964 which were derived 
from radar range and range-rate information on the satellite by the GSFC Data Systems and 
Tracking Directorate. These elements represent a continuous record of "orbit" updates on Syn- 
com n about every week for this four-month drift period. 

Table B-l 


Nineteen Consecutive Sets of Syncom II Orbital Elements, Summer and Fall 1964, 
as Reported by the Goddard Space Flight Center. 






Right 



Epoch, 

Semimajor 

Axis, 

(km) 



Ascension 

Argument 

Mean 

(Yr-Mo-Day-Hr) 

(Universal 

Eccentricity, e 

Inclination, i 
(Degrees) 

of the 
Ascending 

of 

Perigee, cj 

Anomaly, 

M.A. 

Time) 



Node, P 
(Degrees) 

(Degrees) 

(Degrees) 


64-07-04-2.0* 

42202.95 

.00086 

32.541 

312.890 

211.723 

336.434 

64-07-07-3.0 

42203.86 

.00089 

32.537 

312.844 

210.487 

354.275 

64-07-13-17.0 

42203.01 

.00086 

32.502 

312.801 

212.151 

205.972 

64-07-21-21.0 

42202.56 

.00095 

32.516 

312.644 

211.691 

270.784 

64-07-27-16.0 

42202.16 

.00094 

32.502 

312.589 

215.560 

194.968 

64-08-03-17.0 

42201.94 

.00092 

32.476 

312.455 

214.833 

214.528 

64-08-11-1.0 

42201.37 

.00090 

32.446 

312.402 

216.021 

337.265 

64-08-17-19.0 

42201.54 

.00097 

32.409 

312.308 

215.342 

251.656 

64-08-25-10.0 

42201.16 

.00099 

32.443 

312.262 

223.045 

113.049 

64-09-01-10.0 

42201.42 

.00087 

32.397 

312.067 

218.832 

121.245 

64-09-09-14.0 

42200.59 

.00093 

32.373 

311.916 

218.010 

186.557 

64-09-15-12.0 

42201.51 

.00093 

32.362 

311.863 

221.244 

156.571 

64-09-22-10.0 

42201.25 

.00096 

32.302 

311.744 

220,581 

131.010 

64-09-29-6.0 

42202.39 

.00088 

32.327 

311.658 

221.814 

73.456 

64-10-06-5.0 

42202.54 

.00092 

32.320 

311.589 

228.369 

55.663 

64-10-13-0.0 

42203.04 

.00094 

32.312 

311.440 

225.495 

347.171 

64-10-20-16.0 

42203.10 

.00091 

32.289 

311.352 

226.877 

229.785 

64-10-26-16.0 

42204.92 

.00083 

32.251 

311.193 

i 226.984 

232.903 

64-11-02-5.0 

42202.76 

.00090 

32.248 

311.138 

225.015 

73.081 


*This is the first Syncom II epoch after the mean drift rate of the satellite was reduced from -0.7 degree /day to -0.5 degree /day by ground 
commanded, on-board cold gas jet pulsing on 4 July 1964. The set of eighteen updated orbits which follow are complete to 2 November 
1964 to the author’s best knowledge. They have been supplied by the Tracking and Data Systems Directorate of the Goddard Space 
Flight Center through the offices of Dr. Joseph Siry. 
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By using these elements in the particle trajectory program called "Item” (Interplanetary 
Trajectory, Encke Method) at GSFC, elements of semimajor axis and inclination and geographic 
position at the first ascending equator crossing after the epochs of Table B-l have been calculated.* 
The calculated semimajor axes and inclinations at these crossings have been plotted in Figure B-l 
along with a similar set from a comparable single particle trajectory calculated by "Item” start- 
ing from the elements of Syncom II of 64-07-04-2.0 hours U.T. This latter trajectory program 
included longitude earth gravity of amount and orientation close to that determined from the actual 
Syncom n drift data. 

In general it can be said that Figure B-l (and unlisted data on other elements) shows that the 
actual elements of Syncom II have been reasonably well determined throughout the analyzed drift 
period, at least for the purpose of this study. 


♦Listed in Table 1, Syncom II Ascending Equator Crossings and Related Drift-Rate Data, Summer and Fall, 1964, of this publication. 
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r 


32.700 

32.600 - 

©•© 

32.500 - 
32.400 - 

32.300 - 

32.200 

176 

6.6172 

6.6170 - 

- ® 

6.6168 - 

6.6166 - 
5) 

6.6164 - 
6.6162 - 
6.6160 

0 




® 


I 

10 


® . 


© • 


© 


© 


I 


I 


I 


© 


I 


©• 


I 


© 


©< 


© 


© 


I 


♦ © # 


I 


©• • 

© 




_L 




186 196 206 216 226 236 246 256 266 276 286 296 306 316 326 

TIME (da/s) ( from January 0.0, 1964) 


© 


® 


® 


® 




® 


® 


® 


® 


I 

20 


1 

30 


® 


l 

40 


I 

50 


® 


l 

60 


®1 

70 


® 


I 

80 


® 


NOTE: THE UPPER & LOWER TIME 
SCALES ARE IDENTICAL. 


L 


_L 


no 


90 100 

TIME (days) (from July 4, 1964, at2.0hrs.) 
Actual Syncom H Data Estimated by Range and Range Rate Tracking at GSFC 


Orbit Data from a Simulated Sync om H Trajectory Beginning with the Elements of 
Syncom H ( epoch 64-07-04-2.0 hrs . ). Determined by tb’e GSFC " Item " Particle 
Trajectory Program using the Following Earth - Gravity Constants: 


120 


130 


140 


150 


H e =3.98627 x 10 


Ro 

= 6.38388 

X 10 3 

J 20 

=1082.19 

X 

0 

0 

J 22 

=- 1.68 x 

10 6 

J 30 

= -2.285 : 

x icr 6 

J 40 

= -2.212 ; 

x 10' 6 

X22 

= -18.0° 



5 krr)3/ Sec 2 
km /earth radii 


Figure B-l — Orbit inclination and semimajor axis for Syncom II and a simulated 
Syncom II trajectory during the summer and fall of 1964, 
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Appendix C 


Analysis of Simulated Syncom II Trajectories 
(Summer and Fall, 1964) for Sensitivity to the 
Equatorial Ellipticity of the Earth 


Table C-l gives elements and related ascending equator crossing drift data during a simu- 
lated Syncom n trajectory computed numerically (by the "Item” program) beginning with the actual 
elements of Syncom n, epoch 64-07-04-2.0 hours U.T.* The gravity field of the particle trajectory 
(including the fields of the sun and the moon) was exactly that used in the orbit determination pro- 
gram for Syncom n, with the addition of equatorial ellipticity introduced into the geoid 1 " of amount 
and location given by the constants 


J 22 = - 1.68 x 10" 6 


and 


x 22 = - 18° • 

Table C-2 presents ascending equator crossing drift data from a simulated Syncom n trajectory 
calculated with the same parameters as that for Table C-l but without equatorial ellipticity in the 
geoid model. 

The bracketed "mean" longitudes and longitude rates in Tables C-l and C-2 have been cal- 
culated precisely as those in Table 1. A graph of the "observed" (calculated) drift data for the 
longitude gravity trajectory is shown in Figure C-l. * 

For the longitude gravity crossing data, sensitivity to the drift theory expressed by the 
equations 


(k) 2 - - A 22 cos 2(k -\ 22 ) + C x (1) 

or 

(M 2 = c x + c 22 f 1 (X) + s 22 f 2 (\) , (3) 


♦See Table B-l in Appendix B. 

'See Appendix A for a representation of the geoid used in this simulation. 
tSee Figure 1 for trajectory without longitude earth gravity data. 
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Table C-l 


Ascending Equator Crossings and Related Data in a Simulated Syncom II Trajectory* 
With Longitude Dependent Earth Gravity, July-November 1964. 


Crossing 

Number 

<i) 

Date, 

1964 

(mo-day) 

Time.Tj 
(days from 
initial epoch) 

AT, 

T i' T i-l 

(days) 

Longitude,* 1^ 
(degrees) 

AL, 

H-H-i 

(degrees) 

Observed 
Longitude 
Drift Rate,* 
AL/AT 
(degrees/day) 

Inclination 
at Crossing, 
*s 

(32.0°+ data) 

Semimajor 
Axis at 
Crossing, a s 
(6.616 + 
data x 10, ‘ 4 
earth radii) 

Square of 
Observed 
Longitude 
Drift Rate, 

Y (observed) 
(rad/sid day) 2 

Square of 
Computed 
Longitude 
Drift Rate,** 
Y (computed) 
(rad/sid day) 2 

Square of Residual 
Longitude Drift 
Rate, 

Y (obs.) - Y(comp.) 
(rad/sid day) 2 

1 


1.4751 


-171.727 



.540 

4.743 





7-7 



(-172.68) 


(-.4793) 



6.961 x 10 s 

6.903 x 10-5 

5.805 xlO- 7 

2 


5.4693 

3.9942 

-173.642 

-1.915 


.539 

6.850 





7-12 



(-175.05) 


(-.4714) 



6.733 

6.753 

-2.033 xlO' 7 

3 


11.4605 

5.9912 

-176.466 

-2.824 


.519 

5.038 





7-19 



(181.68) 


(-.4644) 



6.533 

6.565 

-3.152x lO' 7 

4 


19.4487 

7.9881 

179.824 

-3.7095 


.512 

5.400 





7-26 



(178.45) 


(-.4582) 



6.361 

6.403 

-4.212x 10' 7 

5 


25.4397 

5.9910 

177.079 

-2.745 


.493 

4.085 





8-1 



(175.72) 


(-.4556) 



6.288 

6.286 

2.234xl0’ e 

6 


31.4306 

5.9909 

174.350 

-2.729 


.485 

3.561 





8-6 



(173.44) 


(-.4550_ 



6.271 

6.202 

6.846 xlO* 7 

7 


35.4246 

3.9940 

172.533 

-1.817 


.473 

5.029 





8-12 



(170.74) 


(-.4484) 



6.092 

6.122 

-3.062x 10~ 7 

8 


43.4123 

7.9877 

168.951 

-3.582 


.452 

2.146 





8-17 



(167.61) 


(-.4467) 



6.045 

6.055 

-9.593 x 10- 8 

9 


49.4032 

5.99085 

166.275 

-2.676 


.439 

4.450 





8-23 



(164.94) 


(-.4445) 



5.986 

6.020 

-3.375x 10* 7 

10 


55.3939 

5.99072 

163.612 

-2.663 


.418 

1.215 





8-30 



(161.83) 


(-.4472) 



6.058 

6.005 

5.3 17 x 10* 7 

11 


63.3817 

7.9878 

160.040 

-3.572 


.399 

3.912 





9-6 



(158.71) 


(-.4445) 



5.987 

6.019 

-3.159x 10' 7 

12 


69.3724 

5.9907 

157.377 

-2.663 


.381 

1.326 





9-12 



(156.04) 


(-.4470) 



6.053 

6.053 

-2.808 x 10" 9 

13 


75.3632 

5.9908 

154.699 

-2.678 


.372 

4.047 





9-19 



(152.91) 


(-.4485) 



6.094 

6.119 

-2.485x 10* 7 

14 


83.3509 

7.9877 

151.117 

-3.583 


.346 

1.373 





9-26 



(149.75) 


(-.4554) 



6.283 

6.214 

6.904 x 10 7 

15 


89.34185 

5.99095 

148.3885 

-2.728 


.336 

4.656 





10-3 



(146.57) 


(-.4564) 



6.312 

6.336 

-2.467 xlO' 7 

16 


97.3297 

7.9879 

144.7425 

-3.646 


.315 

2.466 





10-11 



(142.89) 


(-.4632) 



6.500 

6.510 

-9.880x 10* 8 

17 


105.3178 

7.9881 

141.042 

-3.700 


.301 

5.538 





10-18 



(139.64) 


(-.4692) 



6.668 

6.690 

-2.134 x 10 * 7 

18 


111.3089 

5.9911 

138.232 

-2.811 


.289 

3.514 





10-23 



(137.28) 


(-.4779) 



6.920 

6.834 

8.643xl0- 7 

19 


115.3031 

YT9942 

136.323 

-1.909 


.286 

6.560 





10-28 



(134.89) 


(-.4795) 



6.966 

6.991 

-2.524 x 10 ’ 7 

20 


121.2944 

5.9913 

133.450 

-2.873 


.268 

5.485 





11-7 



(130.02) 


(-.4912) 



7.309 

7.341 

-3.159x10-7 

21 


T3 5.2746 

13.9802 

126.583 

-6.867 


.250 

7.603 





Average ( ' s =32.401°; Average a s = 6.6164238 earth radii 


♦Simulated Trajectory Constants 

^ = 3.98627 x 10 5 km 3 /sec 2 x t2 = -18.0° 

^ = 63.78 3 88 x 10 3 km/earth radii J 30 = -2.285 x 10* 6 
J 20 = 1082. 19 x 10* 6 1 4^ = -2.12x 10" 6 

J 22 = -1.68x 10-6 

The unbracketed longitudes are those printed out from the 
trajectory program. The bracketed longitude rates are 
first difference estimates at the bracketed longitudes. 


♦♦Computed from Equation 3 using the values of C,, C 22 and § 22 
below. Results from the "Least Squares" estimation of the 
parameters C, , C 22 and S 22 in Equation 3 for this simulated 
tr aj ectory: 


C, = 8.4 18 lx 10 ' 5 rad/sid. day 2 S(S 22 ) = 3.328 x 10 
C =-1.4060x10-® 


s 22 

S(C 2J 


1.0392x 10- 6 
) = 3.677xl0- e 


S(standard error of 
estimate) = 4.401 x 10' 7 
rad/sid. day 2 




Table C-2 


Ascending Equator Crossings and Related Data in a Simulated Syncom II Trajectory 
Without Longitude Earth Gravity, Summer and Fall 1964.* 


Crossing 

Number 

(j) 

Time 
(1964) 
(Mo. -Day) 

Time 

(Tj ) 

(Days 

from 

initial 

epoch) 

AT, 

( T j+l“ T i) 

(Days) 

Longitude 

M 

(Degrees) 

AL 

( L i+r L j) 
(Degrees) 

Longitude 
Drift Rate, 
(al/al) 
(Degrees) 
per day) 

i 


1.4751 


-171.737 




7-8 



(-173.2) 


(-.4815) 

2 


7.4665 

5.9914 

-174.622 

-2.885 



7-15 



(-176.5) 


(-.4789) 

3 


15.4549 

7.9884 

-178.448 

-3.826 



7-23 



(179.6) 


(-.4793) 

4 


23.4435 

7.9886 

177.723 

-3.829 



7-31 



(175.8) 


(-.4797) 

5 


31.4319 

7.9884 

173.891 

-3.832 



8-8 



(172.0) 


(-.4802) 

6 


39.4204 

7.9885 

170.055 

-3.836 



8-16 



(168.1) 


(-.4789) 

7 


47.4088 

7.9884 

166.229 

-3.826 



8-24 



(164.3) 


(-.4782) 

8 


55.3973 

7.9885 

162.409 

-3.820 



9-1 



(160.5) 


(-.4812) 

10 


63.3858 

7.9885 

158.565 

-3.844 



9-9 



(156.6) 


(-.4782) 

11 


71.3741 

7.9883 

154.745 

-3.820 



9-17 



(152.8) 


(-.4786) 

12 


79.3626 

7.9885 

150.922 

-3.823 



9-25 



(149.0) 


(-.4801) 

13 


87.3510 

7.9884 

147.087 

-3.835 



10-3 



(145.2) 


(-.4788) 

14 


95.3394 

7.9884 

143.262 

-3.825 



10-11 



(141.3) 


(-.4792) 

15 


103.3278 

7.9884 

139.434 

-3.828 



10-19 



(137.5) 


(-.4787) 

16 


111.3162 

7.9884 

135.610 

-3.824 



10-27 



(133.7) 


(-.4811) 

17 


119.3047 

7.9885 

131.767 

-3.843 



*Initial Orbit Epoch of the Simulated Trajectory (yr.-mo.-day-hr, universal time) 64-07-04-2.0. 
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-7 NOV, 1964 


MEAN DRIFT RATE 
LONGITUDE* 



Goddard Space Flight Centers " Item” Particle Trajectory Input Data 

Initial Elements: Syncom 2 of Epoch 64 - 07 - 04 - 2. 0 hrs 
as Computed at the Goddard Space Flight Center 

Sun & Moon Gravity: Included 

Earth - Gravity Constants: P e - 3. 98627 x 10 5 km 3 / sec 2 

R 0 =6378. 388km 
J 20 = 1082. 19 x 10- 6 
J 22 = - 1.68 x 10- 6 
A 22 = - 18.0° 

J 30 = -2.285 x lO' 6 X 7 JULY, 

J 4 0 — -2.12 x 10' 6 , 1964 


* ACCORDING TO THE 0 
DRIFT RATE THEORY OF 
EQUATION (1): 
COMPATIBLE WITH 
EARTH EQUATORIAL 
ELUPTICIY OF MAGNITUDE: 

J 2 2 = “ 1 . 75x 10” 6 / 

AND ORIENTATION: 

A„ = - 18.2° 


(These are the values used 
at Goddqrd to compute the 
given syncom orbits in 
Appendix A . ) 



6 SEPT, 
1964 


"OBSERVED" CROSSINGS 
FROM THE NUMERICALLY 
INTEGRATED TRAJECTORY 


(Compare with Trajectory Input Data Above) THE LONGITUDE REFERRED TO IS THE LONGITUDE OF THE 

ASCENDING EQUATOR CROSSING OF THE PARTICLE TRAJECTORY. 
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Figure C— ] — Drift rate in a simulated Syncom II trajectory between 
1 87° and 130° East, summer and fall , 1964. 

has been assessed by a revised multiple regression IBM "Share” program* in exactly the same 
manner as the sensitivity of the actual Syncom II data has been tested. Here the "noise sources" 
in the experiment are due just to sun, moon, and earth zonal gravity unaccounted for by second- 
order tesseral earth gravity as well as errors arising from the fact that Equation 3 is an approxi- 
mation of the actual drift perturbation due to such second order longitude gravity. 

The results of this regression analysis (for C 1? C 22 , S 22 in Equation 3) on the longitude gravity 
simulated trajectory are summarized in Table C-l as: 


Cj (sim. ) = 8.4181 X 10“ 


rad 


( s i d day ) 2 

C 22 (sim. ) = - 1.4060 x 10~ 6 

S 2 , (sim. ) = 1.0392 x 10“ 6 


(C-l) 

(C-2) 

(C-3) 


*Kunin, M. J., "Snap-7090 Multiple Regression Analysis Program,” IBM 7090 Program No. 183, "Share” Distribution Number 1289, New 
York: Shell Oil Company Data Processing Department, January 3, 1962. 


28 


DRIFT RATE (-degrees/ day) 


S (^22, sim.) 

3.677 x 1CT 8 

(C-4) 

S (^22,sim.) 

3.328 x 10' 8 

(C-5) 

s £ (standard error of estimate, 

simulated) = 4.401 x 10~ 7 . 

(C-6) 


It is noted that s e (simulated) in this drift period is about a third of s e (actual) during the same 
pei'iod and over about the same number and location of drift rate "observations” * But s e (sim.) in 
the previous Syncom n drift period analyzed (Reference C-l) was only about two-thirds of s e (ac- 
tual) in that period. The greater amount, s e (actual), in both periods was roughly the same. But 
s g (sim.) for the earlier drift period was twice what it is in this new study. It appears likely then 
that the actual "observation" error in this new drift period is somewhat greater in this new four- 
month period than such "position" error in the old. This is the error in the time and crossing 
data supplied by the GSFC Tracking and Data Systems Directorate, as distinguished from sun, 
moon, and other second order gravity noise. Using the values from Equations C-2 through C-4 in 
the equations 


- (C 


2 + o 2 

22 °22 


\ 1/2 


(07) 


and 



(08) 


it is found that 


J 22 ( s im. -calcul at ed) - J 22 (apparent, in the simulated trajectory) D - (1.75 ±0-04) x 10~ 6 (C-9) 

and 

A. 22 ( sim. - calculated) = \ 22 (apparent, in the simulated trajectory) = - (18.2 ±0.6)° . (C-10) 

Since 


J^ 2 (actual, in the simulation) - “1.68 x 10 6 


See Equation 1 3 in this publication. 
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and 


k 22 (actual, in the simulation) - - 18 ° , 

it is apparent that the error biases due to sun, moon and second order gravity "noise" inherent in 
the use of Equation 3 to approximate the longitude gravity drift in this 24-hour satellite arc are 
acceptably small. The good experience with these closely simulated trajectories in References C-l 
and C-2, as well as here, suggests the justification of utilizing the simulated results to predict the 
likely small, accumulated bias in the actual Syncom n data reduction from these most well-known 
or calculable error sources.* Other techniques of removing these calculable errors point for 
point, instead of in this comprehensive manner (for a single drift arc), are being explored. 


REFERENCES 

C-l. Wagner, C. A., "The Equatorial Ellipticity of the Earth from Two Months of Syncom n Drift 
over the Central Pacific," NASA Technical Note D-3315, 1966. 

C-2. Wagner, C. A., "Determination of the Ellipticity of the Earth T s Equator from Observations on 
the Drift of the Syncom n Satellite," NASA Technical Note D-2759, May, 1965. 


•This was done in Equations 20 and 21 in this publication. 
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Appendix D 

List of Symbols 

A 22 Second-order gravity-drift constant associated with the drift regime of a 24-hour satellite. 
C 22 Cosine parameter of earth -equatorial ellipticity. 

C k "Best" or "least -squares" estimation of the parameter c i . 

j 22 Second-order tesseral constants of the earth’s gravity field (associated with magnitude and 
geographic orientation of equatorial ellipticity). 

R 0 Mean equatorial radius of the earth. 

S 22 Sine parameter of earth -equatorial ellipticity. 

a Q Major equatorial radius of the earth. 

a s Semi major axis of the orbit of the near -synchronous satellite. 

b Q Minor equatorial radius of the earth, 
e Orbit eccentricity; when subscripted, equilibrium. 

1 s Inclination of the orbit of the near- synchronous satellite. 
s( ) Estimate of standard deviation of random variable ( ). 
t . Random error associated with observation i . 

i 

k Geographic longitude of the ascending equator crossing. 

<j\ Constant variance of e . (estimate of a e is called s e ) 
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